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Abstract. This lecture reviews fundamental physical processes on star 
formation in galaxy interactions and mergers. Interactions and merg- 
ers often drive intense starbursts, but the link between interstellar gas 
physics, large scale interactions, and active star formation is complex 
and not fully understood yet. We show that two processes can drive 
starbursts: radial inflows of gas can fuel nuclear starbursts, triggered 
gas turbulence and fragmentation can drive more extended starbursts 
in massive star clusters with high fractions of dense gas. Both modes 
are certainly required to account for the observed properties of star- 
bursting mergers. A particular consequence is that star formation scal- 
ing laws are not universal, but vary from quiescent disks to starbursting 
mergers. High-resolution hydrodynamic simulations are used to illus- 
trate the lectures. 



1 Introduction 

This lecture reviews the physical processes of Interstellar Medium (ISM) dynamics 
and Star Formation (SF) in galaxy interactions and mergers. Galaxy collisions 
. . ■ and mergers are one of the fundamental processes of galaxy assembly in A-CDM 
\ cosmology. Mergers are often associated to "starbursts" where the rate of star 
k/, • formation can increase by a larger factor, and the gas consumption timescale 
\ | drops. Here we aim at explaining the observed properties of these merger-driven 
C3 • starbursts with a theoretical point of view and with the help of hydrodynamic 
simulations, and compare to isolated disk galaxies that form stars quiescently. 

The fundamental processes of star formation in any galaxy are firs reviewed in 
a brief and simple way - a thorough review of these is provided by Elmegreen (this 
volume; larXiv:1101.310"8l to 1101.3113). The role of ISM turbulence is highlighted. 
We then study the dynamics of galaxy interactions and mergers, and the general 
properties of SF in mergers. The standard theory for merger-induced starbursts, 
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based on tidally-driven gas inflows triggering a nuclear starbursts, is reviewed 
and its limitations are pointed out - real starbursts tend to be stronger, and more 
extended spatially. We show with modern high-resolution simulations that another 
triggering process is the increase of ISM turbulence and fragmentation, resulting 
from the interactions and inducing an excess of very dense, rapidly star-forming 
molecular gas. Merger-induced starbursts then have a nuclear component, but also 
an extended component in massive and dense star clusters, as generally observed. 
The gas compression is thus both global (over the whole galaxy) from the tidal 
field, and local (on scales of 0.1-1 kpc) from local shocks in the turbulence ISM. 
This explains recent observations, such as a dual law of star formation between 
quiescent disks and starbursting mergers, and an excess of very dense molecular 
gas in Ultra Luminous Infrared galaxies. 

Structure formation and (globular) star cluster formation/evolution, also in- 
cluded in this lecture, a re briefly summari zed at the end of the present text, and 
the reader is referred to Bournaudl ( 2010al ) for more details on these aspects. 



2 Galaxy-scaled Star Formation and ISM dynamics - basic facts 

2.1 Disk stability and star formation 

Star formation on large scales proceeds by gravitational instabilities, which create 
dense clouds in which the gas cools more rapidly, and continue to fragment and 
collapse into start through gravitational (and thermal) instabilities. 

The simplest prescriptions on whether a large (kpc-size) region of a galaxy will 
form stars and at which rate can be obtained by studying its gravitational stability. 
For a gas medium of size L, density p and velocity dispersion^] a, pressure forces 
on the surface vary as L 2 , and necessarily overcome self- gravity, which scales as 
the mass, hence as L 3 . For large-enough scale, gravitational forces eventually 
dominate and large-enough regions are gravitationally unstable (Fig.[T]). 

The critical scale, which is also the most unstable one (the fastest growing 
mode) is the Jeans scale, set by: 



L Jeans * ^ (2.1) 

If gravitational collapse occurs, gas becomes denser, so the cooling time drops 
and the gas effectively gets colder, making pressure forces more and more negligible 
compared to gravity. It is hence realistic to consider the collapse of gas clouds as 
a "free fall" process, i.e. not regulated by pressure forces, at least until the scales 
of magnetic pressure in dense cores are reached, or until star formation has taken 
place and young stars start to heat and ionize the gas. The gravitational free-fall 
time is given by: 



This velocity dispersion can be thermal (microscopic) and/or turbulent (macroscopic). 
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«" = vS (2 - 2) 

At fixed density p in an infinitely large system and without external rotation, a 
large-enough region will always be unstable. But the question to decide on dense 
gas cloud formation and subsequent star formation is whether sufficient mass can 
be gathered over a size not larger than the entire galaxy, and even over a size not 
larger than the disk scale height (above which the mass stops growing as cx pL 3 ). 
Thus a given region of a galaxy may be dense enough or not to be gravitationally 
unstable and form dense gas and stars. 

Another aspect is the stabilization by rotation and shear in galactic disks. Too 
large regions will be torn apart by the shear faster than they tend too collapse, 
namely faster than the gravitational free-fall time. This process stabilizes regions 
tha t are large enou gh to undergo large differential rotation. This is formalized by 
the iToomrel ( 1964 ) parameter: 



Q = (2-3) 

where E is the surface density of the disk, k is the so-called epicyclic frequency, 
without a factor 2 from the angular speed in galactic disks with standard rotation 
curves (3.36 replaces it in the definition of Q for a collisionless fluid, i.e. a stellar 
disk). 

Q > 1 means that any scale is stabilized either by pressure or rotation - scales 
smaller than the Jeans length are stabilized by the gas pressure, and the Jeans 
length is already large enough to be stabilized by rotation. 

Regions with Q < 1 can collapse, fragment, form cold gas clouds and star- 
forming cores therein, and this collapse should take a gravitational free-fall time. 




Fig. 1. An homogeneous region of volume density p is unstable if its size is at least the 
Jeans length, which enables gravity to overcome pressure forces from the internal velocity 
dispersion a. If the region is embedded in a rotation disk, the collapse can occur only if 
the size is small enough compared to the imposed shear ; the Q < 1 Toomre criterium 
indicates that the combination of both constraints can be satisfied. 
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2.2 The heterogeneous and turbulent ISM 

The simple description above could give a simple model to describe star formation 
in galaxies in theory and in simulations. Only regions where Q < 1 can form 
stars, and then a certain fraction of gas is converted into stars within one or a few 
free-fall times, the involved fraction being much smaller than one (typically a few 
percent) because the formation of the first stars will heat the surrounding gas and 
regulate the formation of the next stars (although young stars can also remotely 
trigger the formation of other stars (e.g., Elmegreen fc Lada . 1977?) ). 



However, the involved critical scale in this "Q<1 approach", the Jeans length 
or Toomre length (this two scalelengths are of the same order of value), is of 
the order of 100-1000 pc for normal spiral galaxies such as the Milky Way. If we 
consider a portion of the ISM of a few hundreds of pc, we will find warm atomic gas 
at densities of a few atoms per cubic centimeter and temperatures of a few 10 3 K, 
colder and denser atomic gas, molecular gas with densities of 10 3 to 10 6 cm~ 3 , 
ionized gas at ~ 1 cm" 3 and a few 10 4 K, etc. It is then impossible to define a 
typical density and velocity dispersion (or sound speed) in the ISM on the scales 
that are relevant to measure a Toomre parameter, making this theory impossible 
to apply to real galaxies. 

The only framework in which the basic theory of gravitational instability in 
disks can be directly applied are maybe the resolution-limited SPH simulations 
that are a plenty in the literature since almost two decades. Most of these simu- 
lations (except a few recent ones) do not resolve an heterogeneous ISM, the only 
density fluctuations arising on scales of a kpc and above. There a relevant Toomre 
parameter can be measured and applied as a star formation criterium. But at the 
same time, this simple fact means that the involved ISM model has a quite limited 
realism. 

If we consider the various phases of the ISM in the Galaxy, up to one third of 
the volume may be filled with warm and hot gas above 10 4 K, but at least two 
thirds of the mass (and 100% of the mass involved in dense gas cloud formation 
and star formation) are in denser and colder phases with temperatures between, 
say, 10 K and 10 3-4 K. The sound speed in atomic hydrogen at 10 4 K is about 
10 km s" 1 , thus the star-forming phases of the ISM have a sound speed of the 
order of a km s _1 (0.1 to 10 km s _1 ). 

Observed linewidths of the gas in spiral galaxies are at least a few km s 1 , they 
can reach 50 or 100 km s 1 in interacting galaxies (see details and references later), 
starbursting galaxies, and in gas-rich galaxies at high redshift. This is observed 
even on scales much smaller than the rotational velocity gradient in disk galaxies, 
which shows that supersonic non-circular motions are ubiquitous in the ISM of 
galaxies. 

Interestingly, observations of the density power spectrum of nearby galaxies 
suggest that these random motions follow a cascade process from a large injection 
scale down to dissipat ion on small scales ([Block et all l2010t iDutta et all l2009t 



Elmegreen et all 11993). The spectral signature of t his cascade is just as expected 



from turbulence experiments ([Petersen et all 120071 ). This indicates that the non- 
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Fig. 2. This figure illustrates a density PDF, with a log-normal shape and typi- 
cal value for the galactic ISM. The intrinsic spread of the PDF increases with the 
Mach number of the non-circular (turbulent) motions in the ISM |Wada et all |2002| ; 
iKrumholz fc Thompsonl 120071 ) . The log-normal shape exist is the turbulence is (1) 
isothermal (but in the real non-isothermal ISM the high Mach numbers keep the PDF 
quasi log-normal anyways) in a steady-state (i.e. for an isolated self-regulated disk but 
not necessarily for a suddenly disturbed system in an interaction). 



circular motions in the ISM can be named turbulence. A key property of this 
ISM turbulence is that it is supersonic, since the involved motions largely exceed 
the sound speed. Mach numbers of a few for the cold star-forming ISM phases 
are typical for spiral galaxi es, they may frequently reach ~ 10 in starbursts, and 



high-redshift galaxies (e.g., Forster Schreiber et al. . 20091 ). 



2.3 Star formation in the turbulent ISM 

The connection between ISM turbulence, dense gas phases and star formation 
is best understood by examining the point distribution function of the local gas 
density (density PDF). In its mass- weighted (resp. volume- weighted) version, the 
density PDF represents the mass (volume) fraction of the ISM in bins of density. 
We here use mass- weighted versions. 

In an system at equilibrium (e.g. isolated disk galax y) , supersonic turbulence 
generates a log-normal density PDF (|Wada et all |2002| ) (see Fig. 2). Only small 
mass fractions of the whole ISM are at found number densities below ~ 1 cm~ 3 
or in very high density regions above, say, 10 4 cm -3 , most mass being at ten to a 
few hundreds of atoms per cm 3 (typically observed as HI or CO-traced molecular 
gas). The width of the log- normal PDF depends on various factors, but the main 
dependence is on the turbulent Mach number, with the PDF width a increasing 



when the Mach number increases (see Krumholz fc Thompson! . 12007 . for details) 



A more turbulent ISM will have a larger spread in its density PDF, because the 
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turbulent flows will sometimes diverge and create low-density holes, and some- 
times converge (or even shock) into very dense structures that can further collapse 
gravitationally. Note that log -normal density PDFs are observed, at least at the 



scale of molecular complexes ()Lombardi et al. . 2010h 



The density PDF is a useful tool to describe the star-forming activity of a 
galaxy at a given gas content. For a simplified description, one can consider 
that star formation takes place only in the densest gas phases (i. e. above some 
density threshold), and that in these dense regions the local star formation rate 
follows, for instance, a fixed efficiency pe r fre e-fall time: <9rfr = €ffP^/tff oc 
(see detailed theory in lElmegreenll2002i and lKrumholz fc Thompson] 12007 ). Even 
if the local rate of star formation follows a different behavior than this purely 
density-dependent model (which is physically motivated by the gravitational col- 
lapse timescale), then the first step remains that star formation proceeds only 
in the densest gas phases once dense-enough clouds have formed. The weight of 
this " star- forming phase" is given by the high-density tail of the PDF. Hence, the 
fraction of dense gas along the density PDF is a key parameter for the global star 
formation activity of any given galaxy. The density PDF will thus be a major tool 
to understand starburst activity in galaxy mergers. 



3 Dynamics of galaxy interactions and mergers 

3.1 Dynamical Friction 

When two galaxies start orbiting around each other, the first long-range process 
that drives their evolution is dynamical friction. When one object (e.g., the stellar 
body of a galaxy) flies through another mass distribution (e.g., the dark matter 
halo of the other galaxy), the gravitational forces will unavoidably create a mass 
accumulation from the latter system along the wake of the first object. This will 
in turn induce braking forces. This process also acts on non-overlapping systems, 
although it is much weaker then. Dynamical friction is what drives the merger 
of interacting galaxies into single, more massive objects. In the process a small 
fraction of the mass is expulsed and becomes unbound, carrying away a substantial 
fraction of the initial energy at high velocity, so that the merger remnant is a lower- 
energy, usually more compact system. 



3.2 Violent Relaxation 

The typical merger timescale for a pair of massive galaxies is about a billion year, 
although large variations are possible for small impact parameters or high en- 
counter velocities. This is barely larger than the rotation period of the outer disk of 
a spiral galaxy. Thus, stars (or gas clouds or dark matter particles) in an interact- 
ing and merging system undergo a variation in the mass distribution/gravitational 
potential that is relatively rapid compared to their own dynamical timescale. In 
such conditions stellar orbits do not follow a slow secular evolution anymore, but 
can be completely redistributed. For a given star, the average energy along its 
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final post-merger orbit can be quite different from the initial average energy along 
the initial pre-merger orbit: this large difference arise because of the rapid change 
in the mass/potential distribution and defines the violent relaxation. Some stars 
(or gas clouds or dark matter particles) will then undergo large decrease of their 
average energy, which implies that they sink towards the central regions: the cen- 
tral regions of galaxies thus become denser and more compact in a merger. Other 
stars will gain energy and increase their orbit semi-major axis - some even escape: 
a more extended outer envelope in the mass distribution results. 

When one writes the projected radial surface density profile of a galaxy as a 
Sersic profile log(S) ~ fj, ~ e~(r — a) 1 /"/^), where n is the Sersic index, spiral 
galaxies have exponential disks with n ~ f . A merger remnant will have a steeper 
central profile and a more extended outer envelope, which corresponds to n > 1. 
The high Sersic indices of elliptical galaxies (say, 2 < n < 6) are thus consistent 
with the common idea that they are the end-product of spiral-spiral mergers. 



3.3 Gas response and tidal torques 

The tidal field A first driver of the gas response in galaxy interaction is the 
tidal field. Typically, the tidal field of a single galaxy is disruptive at large radius 
and outside its own disk, but can be compressive in the inner regions if the mass 
distribution has a low Sersic index. During a distant interaction or the early stages 
of a merger, a given galaxy is then mostly affected by a disruptive tidal field from 
its distant companion. More advanced mergers can have c ompressive tidal field 
over large regions, when the galaxies begin to overlap (e.g., Renaud et al. . 20081 
for the Antennae). 

While a disruptive tidal field naturally tends to expel material from disk galax- 
ies, it is not the only mechanism responsible for the formation of long tidal tails, 
and it cannot explain central gas inflows: this is in fact mostly driven by gravity 
torques. 



Gravity torques, inflows and outflows The tidal field from a companion 
breaks the symmetry of the gravitational potential. This induces a response of 
the disk material, in particular its cold gas, which is more pronounced for pro- 
grade! interactions. The gas can form an interaction-driven pair of grand - design 
spiral arms, like modeled in detail for instance for M51 bv iDobbs et al. ( 2010f ) 



- note though that grand-desig ns spiral do not require an interaction to form 
(lElmegreen fc Thomasson, 19931 ). and are not connected to interactions in obser- 
vations ( van den Berghl . 20021 ). The gas response is often more complex than a pair 
of spiral arms, but a general feature remains: inside the corotation radiu^l, the gas 
preferentially concentrates on the leading side of the valley of potential; outside of 
the corotation it concentrates on the trailing side. A similar but better-known gas 



2 in which the orbit and disk rotations have similar orientations 

3 The corotation radius is the radius at which the rotation speed of the studied disk equals 
the orbital speed of the companion. 
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Fig. 3. Simplified description of the gas response in a galaxy interaction. Inside the 
corotation radius, gas concentrates on the leading side of the valley of potential and 
undergoes negative gravity torques, which drives a central gas inflow. At the opposite, 
positive gravity torques dominate outside the corotation, which drives lot of the outer 
disk material in the so-called "tidal" tails. In detail the gas response is more complex 
because the system evolves rapidly in a non-equilibrium configuration, but the main 
response is the one shown here. The same response occurs in barred galaxies in a weaker 
but more directly observable way (see Figure [4]) . 



Fig. 4. In a barred galaxy, the valley of potential along the bar induces the same gas 
response as in an interaction (but at a weaker degree). The gaseous arms on the leading 
side of the bar inside the corotation radius, attested by two sharp dust lanes on this 
optical image, are more visible as the system evolves in a steady state. 



response is found in barred galaxi es, where the disk symmetry is broken about the 
same way by the barred pattern (jCombes fc Gerinl . Il985t lAthanassoulal Il992f ) . A 
corotation does not always exist and can move with time, but is in general located 
at a few kpc of radius in the disk. The subsequent process is illustrated on Fig- 
ure [3] inside the corotation, the gas undergoes negative gravity torques, and loose 
angular momentum in a rapid central gas inflow towards the central kpc or less. 
In the outer disk, gas would gain angular momentum and fly out to larger radii in 
long tidal tails. 
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It is worth noting that: 

• The so-called "tidal" tails do not just result directly from the tidal field. 
Gravity torques induced by the interaction are an important process in tidal 
tail formation. If the tidal field was the main factor, few mergers would actu- 
ally have a central gas inflow and related central starburs t. This was already 
illustr ated by the restricted three-body simulations by iToomre fe Toomra 
(Il977h : the tidal field was included but not the full disk response, tidal tails 
were forming but relatively short and no central gas inflow was induced. 
The role of disk self-gr avity and grav itational torques in amplifying the disk 
response was noted in iToomrel (|l98ll ) 



• Stars in a disk galaxy behave in a roughly similar way, but their much larger 
velocity dispersion and collisionless dynamics makes the gravity torques 
much weaker. This is why stars pre-existing to a galaxy interaction are 
generally not found in great amounts in tidal tails (e.g. Due et al. 1997, but 
see for instance Due et al. 2000 in Arp245). For the same reason, there is 
no significant central inflow of pre-existing stars (while the inflow of gas can 
form a lot of new stars in the center). 

Collisional ring galaxies form through a different dynamical process. The central 
impact of a companion onto the disk triggers large epicyclic motions around the 
initial orbits of stars and gas clouds. The superposition of motions and the varia- 
tion of the epicyclic frequency with radiu s result in an expanding way e of ma terial. 
This process is thoroughly reviewed bv lAppleton fc Struck-Marcelll (j 1996h . The 
ring wave is initially an unbound feature, like the spiral waves in an isolated disk, 
but the ring c an become self-gravitat ing, bound, when enough mass is gathered in 



the ring (e.g.. iBournaud et all 120071 ). While the initial process differs, collisional 
rings and are relatively similar to tidal tails in terms of substructure formation, 
so they will hereafter be treated as the same kind of "tidal" structures. 

4 Star formation triggering in interactions - From mergers to star- 
bursts 

4.1 Main observed facts 

Observations of a starburst merger connection The main and simplest 
observation probing the link between starbursts and mergers is that almost all 
starbursting galaxies in the nearby Universe are mergers. In particular, (Ultra- 
luminous Infrared Galaxies (LIRGs and ULIRGs) at l ow redshift are almost 
exclusively strongly interacting and merging systems (e.g., Due et al. . 19971 ). Note 



that at high redshift (z ~ 1 and above), LIRGs are not necessarily mergers because 
"normal" iso l ated d i sk galaxies are gas-r ich enough to reach LIRG-like activity 
( Dadd i et al. . l2010at rTacconi et all 120101) . but it remains true that the strongest 



starbursts (ULIRGs and Hyper-LIRGs) are predominantly merging systems (e.g., 



lElbaz fc Cesarskvl 120031 ) 
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It has been long known that mergers can trigger bursts of star formation 



( Sanders et all 119881 ). But the fact that the most actively star forming systems at 
any given redshift are mostly mergers (LIRGs at z ~ 0, ULIRGs at z ~ 1) does not 
necessarily imply that mergers are always starbursting, or that all mergers even 
undergo a significant starburst phase - we here of course consider "wet" mergers 
where at least one galaxy contains a substantial amount of gas, not "dry" mergers 
of gas-free galaxies. Understanding whether all mergers are starbursting requires 
to be able to select mergers (at least major ones) in large surveys including star 
formation tracers. The main concern i s that any techn i que selecting mergers ba sed 



on morphology (for instance those by IConselice et al.l . l2003t lLotz et all 120081 ) re- 



mains ambiguous, in particular at hig h redshift where morpho logical disturbances 



and asymmetries can arise internally ( Bournaud et al. . 20081 ). Surveys including 



spatially-resolved gas kinematics are certainly a m ore robust tool to id entify on- 
going major mergers and (relatively) isolated disks ( Shapiro et all I2008T ) but such 
samples remain too limited to study the me rger-starbu r st connection. Neverthe- 



less, the most recen t and detailed studies ( Bell et al. . 2006t Jogee et al. . 20091 



Robaina et al. , 2009) are broadly in agreement that mergers do trigger star for- 
mation activity, and specific star formation rate (i.e., SFR divided by the stellar 
mass already in place) being increased by factor of about 3-4, on average: strongly 
starbursting mergers do exist but most mergers do not increase the SFR by a 
factor above about 10 at the peak, likely because strong starbursts are rare, or 
they are short phases compared to the duration of an interaction/merger until 
final relaxation of the post-merger early-type galaxy. 

The factor 3-4 of increase in the SF activity is the average observed value for 
randomly interactions/mergers at any instant of the merging process, not at the 
peak of a starburst activity, so it likely indicates that peak SFR values are higher 
than this. Also, it is the typical measure for "major" interaction/mergers, includ- 
ing mass ratios of, say, 1:1 to about 3:1. There is no evidence of a variation with 
redshift, suggesting that the typical SF activity of "isolated" disks and merging 
systems increase with redshift in roughly similar proportions. 



Morphology of merger-induced star formation Merger-induced star for- 
mation consists, for a large part, of nuclear starbursts taking place in the central 
100-1000 pc (diameter). Nuclear activity in merger-induced st arbursts ha s long 



been emphasi z ed in studies of the merger-SF connection (e.g.. iKeel et al 



Sanders et~ail Il988t iDuc et all Il997t ISoifer et ail . 1 1984 lLawrence et al 



1985; 
1989). 



Nevertheless, the importance of nuclear starbursts among merger-induced star for- 
mation in general has long been overstated in the literature. This may be partly 
because theorists have long been able to explain only nuclear starbursts, or be- 
cause of using infrared and/or Ha luminosities to trace star formation - these are 
dependent on metallicity, dust properties, and threshold effect making them more 
sensitive to compact central starbursts than to spatially extended star formation 
with lower surface density. 

Actually, there is a large fraction of systems in which merger-induced star for- 
mation is spatially extended, taking place outside the central kpc, and in these 
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system the extended star forming component is not just a relics from extended 
SF in pre-merger spiral galaxies, but participates to the total starburst activity, 
even if there is also intense SF in the nucleus. The relatively large spatial extent 
of SF and the problem it poses f or theoretical models (as we will review below) 
was probably first recogniz ed by iBarnesl 12004 in the Mice (NGC 4676A+B, see 
also Chien fc Barnesl I2OI0I in NGC7252). A well-know example if actually the 
Antennae system (NGC 4038+4039), where a large fraction of the starburst ac- 
tivity takes place at several kpc from the nuclei in big "super star cluster" (SSCs) 
formed along density waves throughout t he interacting disk s, and in the shocked 



overlap region between the two gas disks (|Wang et all . 120041) - these star-forming 



components are more important than the dust-enshrouded nuclear star formation 
revealed in the infrared. There are many other examples of interacting/merging 
systems with spatially extended starbursts whose total SFR is far from being domi- 
nated by the nuclear SF acti vity, but rather by ex tended star-forming com ponents: 
Apr 140 (|Cullen et al.l . l2006l ); NGC2207+IC2163 (|Elmegreen et all l2006h in which 
the starburst is almost exclusively non-nuclear with an SFR of 15 M Q yr -1 for 
global gas mass/size about similar to typical Milky Way- type spirals (i.e., a sig- 
nificant increase of SFR likely results from the interaction, but only in a spatially 
extended mode); the Antennae: the C artwheel system - see also many tidally in- 
teracting pairs in lSmith et al.l ( 2008 ). Hancock et al. ( 20091 ). etc. There is to my 
knowledge no thorough census of the contributions of nuclear and extended star 
formation to merger- induced activity, probably because unveiling dust-enshrouded 
nuclear SF and detecting spatially extended SF at lower surface densities requires 
to employ different techniques at different wavelength. Nevertheless, eye examina- 
tion of systems with UV imaging available in samples of interactions and mergers 
(e.g., the Arp Atlas) suggest that spatially extended (>lkpc) SF and nuclear ac- 
tivity (within the central kpc) each contribute to a rough 50% of star formation 
in mergers. 

ISM properties in interacting starbursts Two fundamental properties of 
the interstellar gas in star-forming mergers majors will hold a fundamental role in 
our theoretical interpretation of the starburst activity: 



• gas velocity dispersions are higher in interacting/merging systems than in 
isolated galaxies. At redshift zero, the cold ISM of spiral galaxies has a 
typical velocity dispersion of, say, 5-10 km s , and only gas-rich disk at 
high redshift have higher dispersions. Major mergers can have g as velocity 
dispersions of several t ens of km s _1 , as firs t observed by Irwinl (0*9 94) and 



Elmegreen et al.1 (jl995l ) (see also examples in lBournaud et al.l l2004 although 



limited to star- forming clumps in Ha). These increased dispersions do not 
result from blending of several components in a large telescope beam (a.k.a. 
beam smearing), as they are also observed at high resolution, but trace a 
physical velocity dispersion on small scales. As this applies primarily to gas 
components observed through HI or CO lines, the thermal dispersion (sound 
speed) is lower than 10 km s _1 (for temperatures below 10 4 K), so most of 
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the observed dispersion in isolated galaxies and merging systems consists in 
supersonic gas tur bulence, which usually dominates the cold star forming 
phases of the ISM (|Burkertl 1200(1 (a warm, thermally-supported phase can 



fill a large volume fraction, but involves only a minor mass fraction of the 
ISM). Hence, the (supersonic) ISM turbulence is typically increased by a 
factor of a few in galaxy interactions and mergers. We will discuss later, 
using numerical and theoretical models, that this is probably a trigger of 
increased star formation (through compressive regions in the turbulent flow) 
rather than a consequence of the star forming activity itself (through stellar 
feedback impacting the surrounding gas). 

• Various density phases of the ISM are not filled and distributed in the same 
way as in isolated disk galaxies. There is often a spatial segregation between 
the moderate density, atomic phase, mostly found in the outer regions, and 
the dense r molecular phas es concentrated towards the few central kilopar- 



secs (e.g-. lDuc et all 11997( 1 . More systematically, starbursting systems in the 
nearby Universe with ULIRG-like infrared luminosities (which are almost ex- 
clusively merger-i nduced starbu r sts, se e above) exhibit enhanced HCN/CO 
luminosity ratios ((Juneau et al. , 2009). HCN is usually considered to be a 
tracer of the very dense gas (say, > 10 5-6 cm" 3 ) directly involved in actual 
star formation, while CO lines also (and mostly) trace lower-density regions 
(down to, say, 100 cm~ 3 ) that belong to star- forming clouds but are not 
directly star forming: hence observations point towards an excess of high- 
density gas in starbursting mergers compared to the amount of moderate- 
density molecular gas traced by CO. 

As detailed in the next Section, if the density PDF of the ISM is driven by the 
supersonic turbulent support, then the increased velocity dispersion observed in 
mergers can directly result in a larger spread of the density PDF, which can result 
in an excess of high-density molecules compared to the total molecular gas mass. 
We will illustrate later, using hydrodynamic simulations, how these processes act 
in concert to trigger the star formation activity of interacting and merging galaxies. 

4.2 Standard theory of merger-induced starbursts: gas inflows and nuclear 
starbursts 

In an interacting galaxy pair, the gas distribution becomes non-axisymmctric, 
which results in gravitational torquing of the gas as detailed previously. Gas 
initially inside the corotation radius (typically a radius of a few kpc) undergoes 
negative gravity torques and flows inwards in a more and more concentrated central 
component (inside the central kpc) . Any model for star formation will then predict 
an increase of the star formation rate (global Schmidt-Kennicutt law, models based 
on cloud-cloud collisions, etc). The result is thus a centralized, nuclear starburst. 
As the driving process is gravitatio nal torquing, early restrict ed three-body models 
could already describe the effect ( Toomre fc Toonirel 1972 1. Later models have 



added extra physical ingredients leading to more accurate predictions on the star 
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formation activity: self-gravity (Barnes fc Hcrnquist, 1991), hydrodynamics and 



feedback processes (jMihos fe Hernquis tl. ll996t ICoxet all l2008h . etc. 



A large library of SPH simulations of galaxy mergers, in which the driving 
process is mostly the one presented above (tidal torqui ng of gaseous disk) was 
performed and analyzed by (|Di Matteo et alll2007ll2008l ). This study highlighted 
various statistical properties of merger- induced starbursts. In particular, it showed 
that some specific cases can lead to very strong starbursts with star formation 
rate (SFRs) are increased by factors of 10-100 or more, but that on average the 
enhancement of the SFR in a random galaxy collision is only a factor of a few (3-4 
being the median factor) and only lasts 200-400 Myr. These results were confirmed 
with code compari sons, and found to be in dependent of the adopted sub-grid model 
for star formation ( Di Matteo et al ., 2008). Models including an external tidal field 
to simulate the effect of a large group or cluster found that the merger-induced 
starbursts could be somewhat more efficient in s uch context - but the SFR i ncrease 
remains in general below a factor of 10 or less (jMartig fc BournaudL 120081 ) . 



4.3 Limitations of the standard theory 



The intensity of merger-induced starbursts Numerical simulations repro- 
ducing the interaction-induced inflow of gas and resulting nuclear starbursts can 
sometimes trigger very strong starbursts, but in general the SFR enhancement 
peaks at 3-4 times the sum of the SFRs of the two pre-merger galaxies. This 
factor of 3-4 seems in good agreement with t he most recent observational esti- 
mates ( Jogee et al. , 20091 Robaina et al. . 20091 e.g.). In fact, there is substantial 
disagreement: the factor 3-4 in simulation samples is the peak amplification of 
the SFR in equal-mass mergers. In observations, it is the average factor found at 
random (observed) instant of interactions, and in mergers that are "major" ones 
but not strictly equal-mass ones. Given that typical duration of a merger is at 
least twice longer than the starburst activity in th e models, and th at unequal- 



mass mergers make substantially weaker starbursts (|Cox et all 120081 ). one would 



need a peak SFR enhancement factor of about 10-15 (as measured in simulations) 
to match the average enhancement of 3-4 found in observations. There is thus 
a substantial mismatch between the starbursting activity predicted by existing 
samples of galaxy mergers, and that observed in the real Universe - although the 
observational estimates remain debated and may depend on redshift. 

The spatial extent of merger- induced starbursts The usual explanation of 
merger-induced starbursts accounts only for nuclear starbursts. While these are 
relatively frequent among mergers, they do not necessarily dominate interaction- 
triggered star formation, and spatially extended are actually common (see above, 
and two examples of Figure [5]). Quantitative comparisons of the extent of star 
formation in observations to that pred ic ted by "standard" mo dels have shown 
a significant disagreement iBarnesI (|2004[ ): IChien fc Barnes! (|2010l) . These authors 
highlighted the failure of standard models to account for the extended morphology 
of (at least some) merger-induced starbursts. They also suggested that a sub-grid 
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V 




Fig. 5. Two examples of starbursting interactions/mergers (the Antennae and the 
NGC2207+IC2163 pair), where a large part of the starburst activity is not limited to 
the central kpc regions, but also takes places in super star clusters and giant HII regions 
throughout the disturbed disks, and in overlap regions between the interacting galaxies. 



model of shock-induced star format ion may better acc ount for the spatial extent 
of merger- induced activity (see also Saitoh et al. 20091) . 



4.4 Star formation in mergers: effect of the turbulent, cloudy ISM 

The usually invoked mechanism of gas inflows for merger-induced star formation, 
as reproduced in low-resolution simulations, certainly takes place in real mergers - 
signs of nuclear merger-driven starbursts are a plenty. But it seems impossible to 
explain the typical intensity of merger-induced starbursts and their often relatively 
larger spatial extend, based on these "standard" models. 

Traditional SPH simulations model a relatively warm gas for the ISM, be- 
cause the limited spatial resolution translates into a minimal temperature under 
which gas cooling sh ould not be modeled (it would generate artificial instabilities, 
ITruelove et al. 1997 ). Modeling gas cooling substantially below 10 4 K requires 
"hydrodynamic resolutions'^ better than 100 pc. Cooling down to 100 K and 
below can be modeled only at resolutions of a few pc. The vast majority of 
existing merger simulations hence have a sound speed of at least 10 km/s and 
cannot explicitly treat the supersoni c turbulence that characterizes most of the 
mass in the real ISM (jBurkerti 120061 ) . Turbulent speeds in nearby disk galaxies 
are of 5-10 km/s, for sound speed of the order of 1-2 km/s in molecular clouds, 
i.e. turbulent Mach numbers up to a few. These are even higher in high-redshift 



4 i.e., the average smoothing length in SPH codes, or the smallest cell size in AMR codes 
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disks (For ster Schreiber et al. , 2009), and in mergers (references above), but not 



necessarily for the same reason. 

Increased ISM turbulence in galaxy mergers is thus absent from the model- 
ing used in most hydrodynamic simulations to date. Some particle-based models 
have nevertheless been success f ul in reproducing these incr eased gas dispersions 



( Elmegreen et all Il993t IStrucld . 119971 : iBournaud et all 120081 1 , indicating that it is 
a consequence of the tidal interaction which triggers non-circular motions, rather 
than a consequence of starbursts and feedback. It should then arise spontaneously 
in hydrodynamic models, if these a capable of modeling gas below 10 3_4 K. 



High-resolution simulations with multiphase ISM dynamics Adaptive 
Mesh Refinement (AMR) codes allow hydrodynamic calculations to be performed 
at very high resolution on adaptive-resolution grids. The resolution is not high 
everywhere, but the general philosophy is to keep the Jeans length permanently 
resolved until the smallest cell size is reached. That is, the critical process in 
the collapse of dense star-forming clouds, namely the Jeans instability (or Toomrc 
instability in a rotating disk) is constantly resolved up to a typical scale given by 
the smallest cell size, or a small multiple of it. AMR simulatio ns of whole galaxies 



2009L 



have recently reached resolut ions of a few pc for disk galaxies (jAgertz et al. 
e.g.,), and even 0.8 pc lastly (jBournaud et all l2010cj ). Such techniques have been 



first employed to model ISM dynamics and star formation in galaxy mergers by 
Kim et al. | (<2009h : iTevssier et al.1 (|2010l) . The main physical improvement is not 
necessarily directly the higher spatial resolution, but the possibility to model gas 
cooling down to low temperatures (below 300 K) while still resolving the typical 
Jeans length of gravitational instabilities. In theory, SPH techniques could also 
allow low temperature floor at high spatial resolution, but in practice this may 
be more expansive if not prohibitive, a nd has not been done for mergers and star 
formation (but see lHopkins et al. 2011 but high- resolution, low-temperature floor 
SPH disk models). 

We here illustrate the results of such simulations using a sample of AMR sim- 
ulations of 1:1 mergers of Milky Way-type spiral galaxies, performed with a res- 
olution of 4.5 pc and a barotropic cooling model down to ^50 K, tech nically 



similar to the isolated disk simulation described in ( Bournaud et al. . l2010d) . Star 



formation takes place above a fixed density threshold and is modeled with a lo- 
cal Schmidt law, i.e. the local star formation rate density in each grid cell is 
Psfr = e ffP/tff P 1 ' 5 where p is the local gas density and tff the gravitational 
free-fall time, and the efficiency e// is a fixed parameter. Supernova feedback is 
included. Further details and results for whole sample will be presented in Powell 
et al. (in preparation). An individual model of this type (but at slightly lower 
resolution and without feedback), matching the mor phology and star for mation 
properties of the Antennae galaxies, was prese nted in Teyssier et al. 2010l . More 
details on the simulations can also be found in (|Bournaudil2010bl) . 



In the models, the pre-merger isolated spiral galaxies spontaneously develop 
ISM turbulence at a about 10 km s^ 1 under the effect of gravitational instabilities 
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(and/or feedback), and most star formation takes place in dense complexes of 
dense gas along spiral arms. In some sense, the large-scale star formation process 
is not entirely sub-grid anymore in these simulations, as the first steps of star 
formation, namely are the development of ISM turbulence and the formation of 
dense molecular gas clouds in this turbulent ISM, are explicitly captured - the 
subsequent steps of star formation at smaller scales, inside the densest parts of 
these cold clouds, remain sub-grid. 

A merger simulation is shown in Figure [5] The mass- weighted average of the 
gas velocity dispersion reaches 30 — 40 km/s (Fig. [5]). This strong turbulence 
is consistent with the observations reviewed above. It induces numerous local 
shocks that increase the local gas density, which in turn triggers the collapse of 
gas into cold clouds. Also, gas clouds become more massive and denser than in 
the pre-merger spiral galaxies. The fraction of gas that is dense-enough and cold- 
enough to form star increases, and the timescale for star formation in these dense 
gas entities (the gravitational free-fall time) becomes shorter. As a result, the 
total SFR becomes several times higher than it was in the pre-interaction pair of 
galaxies. The standard process of merger-induced gas inflow towards the central 
kpc or so is also present, but the timescale is substantially longer, so this process 
dominates the triggering the star formation by enhancing the global gas density 
only in the late stages of the merger. 

An example of star formation distribution is shown in Figure |SJ Two conse- 
quences of modeling a cold turbulent ISM in merging galaxies are: (1) the peak 
intensity of the starburst can become stronger (as shown by Teyssier et al. 2010 
for a model of the Antennae) although this is not a systematic effect, and (2) 
the spatial extent of star formation in the starbursting phase is larger. The ra- 
dius containing 50% of the star formation rate (half-SFR radius) can more than 
double. This is because increased ISM turbulence is present throughout the sys- 
tems (see Fig. [5J): this triggers, through locally convergent flows and shocks, the 
formation/collapse of efficiently star- forming clouds even at relatively large dis- 
tances from the nuclei. At least quantitatively, these results put the models in 
better agreement with observations. In these models, the increased ISM turbu- 
lence is also obtained without feedback, showing that it is not a consequence of 
the starbursting activity, but is rather driven by the tidal forces in the interaction. 

4.5 From triggered turbulence to dense gas excess and starbursts 

Modern models are thus capable of resolving ISM turbulence and SF clouds in disk 
galaxies and starbursting mergers. ISM turbulence is self-regulated in isolated disk 
galaxies, and is significantly triggered in merging systems. 

Our simulated of disk galaxies have quasi log-normal PDFs (see example on 
Fig. [7]). The spatial resolution limit converts into a density limit at which the 
density PDF is truncated, and which corresponds to the smallest /densest entities 
that can be resolved. The maximal density resolution in the models presented 
here is around 10 6 cm" 3 , but parsec-scale resolution could capture even higher 
densities (as in Bournaud et al. 2010). 
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Fig. 6. Radial profiles of the gas velocity dispersion (turbulent speed) in an isolated 
galaxy simulation, and three major mergers between two disks of the same initial type. 




Fig. 7. Density PDFs of an isolated disk (left) and a starbursting major merger (right), 
in AMR simulations. The excess of very dense, star-forming gas results from the non 
log-normal shape of the PDF, induced by the rapidly increased turbulence and numerous 
local shocks compressing the gas, more than by a global shift of the initial PDF. 



Representative density PDFs are shown for merger models on Figure[7|(see also 
Fig. \§§ for a moderately starbursting merger (SFR increased by a factor ^3 — 4 
compared to the two pre-merger galaxies taken separately), and a stronger merger- 
induced starburst (factor ~ 10). The PDFs of merging galaxies have a larger half- 
maximum width than those isolated disk galaxies, as expected from the higher 
turbulent speeds that result from the tidal interaction. Also, these PDFs are not 
log-normal, as would be the case for isolated disks with a self-regulated turbulence 
cascade in a steady state. Rapid perturbations of the velocity field in mergers 
result in a substantial excess of dense gas in the density PDF. Such density PDFs 
naturally imply high SFRs, independently from the local SFR prescription, since 
the fraction of efficiently star-forming gas is high. 



High fractions of dense gas in mergers were already proposed by ((Juneau et al 



2009), based on detailed post-processing of merger simulations ai med at re-construc ting 



dense molecular gas phases not resolved in these simulations (jNaravanan et al 



2010h . Here we obtain a qualitatively similar conclusion using simulations that 



18 



Title : will be set by the publisher 




Fig. 8. View of the very young stars (proxy for the instantaneous star formation rate) 
in a major merger model from Teyssier et al. (2010). 



explicitly resolve turbulent motions, local shocks and small-scale instabilities in 
cold ISM phases. This excess of dense gas should have signatures in molecular 
line ratios. If, in a rough approach, we assume that low-J CO lines are excited 
for densities of 100 cm -3 and above, and HCN lines for densities of 10 4 cm~ 3 
and above, then the HCN/CO line ratios could be ~ 5-20 times higher in the star- 
bursting phases of major mergers. Simulations with an somewhat higher resolution 
would actually be desirable to accurately quantify the emission of dense molecular 
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Fig. 9. Statistical distribution of the SFR/HI ratio (measured as the ratio of gas denser 
than 10 5 cm -3 as a proxy for SFR and gas less dense than 10 2 cm -3 as a proxy for HI (as 
illustrated on fig. [7J), for isolated spirals and mergers at different stages. This quantified 
the excess of dense gas illustrated earlier by density PDFs. 

tracers, and opacity effects may also affect the observed HCN/CO ratios. 
4.6 Implications for star formation "laws" 

The interpretation of merger- induced starbursts proposed from our high- resolution 
models is that it is not just a global gas inflow that increases the average gas density 
and increases the SFR, but also that there are strong non-circular motions, high 
turbulent velocity dispersions, causing many small-scale convergent flows and local 
shocks, that in turn initiate the collapse of dense star-forming clouds with high 
Jeans masses. The former "standard" process does take place, but the later can 
be equally important especially in the early phases of mergers. 

We here note £ gQfj the average gas surface density of a galaxy. This is the quan- 
tity that observers would typically derive from the total gas mass and half-light 
radius, or similar quantities. The second mechanism above is a way to increase the 
SFR of a system, and its SFR surface density T,sf r, without necessarily increas- 
ing its average S gas . Actually in our merger models T, gas does increase (as there 
are global merger- induced gas inflows), but £sf.r increases in larger proportions 
(as the starburst is not just from the global merger-induced inflow but also from 
the exacerbated fragmentation of high-dispersion gas) . Going back to the density 
PDFs shown previously, one can note that the fraction of very dense gas (say, in 
the ~ 10 4 ~ 6 cm" 3 range) can increase by a factor of 10-20 in mergers while the 
average surface density E gas increases by a factor 3-5 (see also on Figure 4). As 
a consequence, the £sf.r activity of these systems is unexpectedly high compared 
to their average surface density E gas . 

Figure 4 shows the evolution of a system throughout a merger simulation in 
the (£ gas ;Eg^i?) plane. While our pre-merger spiral galaxy models lie on the 
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standard Kennicutt relation, starbursting mergers have high Ss_f i?/S gQS ratios. 
This is in agreement with observational suggestions that quiescent disks and star- 
bursting mergers do not follow the same scaling relations fo r star format i on, but 
could actually disp lay two different star formation "laws" (jDaddi et all [2010b; 
iGenzel et all l2010h . The offset between the disk and merger sequences proposed 
by Daddi et al. (2010b) is quantitatively recovered in our simulations (Figure 4). 
Post-starburst, post-merger systems lie back on the quiescent sequence, or even 
somewhat below it: these systems contain some dense gas which is somewhat 
stabilized by the stellar spheroid. This is another example of a "morphological 
quenching" effect in early- type galaxies ( Martig et al. . 2009() . and the location of 
our post-merger early-type galaxies in the (E gas ; Ssf_r) diagram may be consis 



tent with observations of nearby ellipticals ([Crocker et all 1201 It ISaintonge et al 



20111 ). without having a substantial effect o n the formation of the red sequence 



at least at low redshift (jFabello et all 120111 ) 



The proposal that disks and mergers follow two different regimes of star for- 
mation by Daddi et al. (2010b) and Genzel et al. (2010) relies for a part (but 
not entirely) on the assumption that different CO luminosity-to-molecular gas 
mass conversion factors apply in quiescent disks and starbursting mergers. In- 
terestingly, our simulations recover the two regimes of star formation without any 
assumption on such conversion factors since gas masses are directly known. But at 
the same time, excess of dense gas found in these merger models suggests that the 
excitation of CO lines would naturally be higher in mergers/starbursting phases 
(although this needs to be quantified in the models), which would mean that the 
assumption of different conversion factors by Daddi et al. and Genzel et al. could 
be physically justified. High molecular gas excitation in Sub-Millimeter Galaxies 
(SMGs, Tacconi et al. 1998) could then naturally result if these are starbursting 
major mergers w i th high gas surface dens ities and a clumpy turbulent ISM (e.g. 
iNarayanan et al. . 2010 : iBournati d. 2010b). However models do not predict a bi- 
modality. Two regimes of star formation appear when mergers are only selected at 
the peak of their starburst activity, but merging systems spend over half of their 
time below the " starburst regime" on the diagram shown in Figure 1101 



4. 7 Summary 

The triggering of star formation in merger largely results from tidal inflows driven 
by the interaction, which fuels a relatively concentrated or "nuclear" starburst. 
This is the main mechanism highlighted by resolution-limited simulations for two 
decades. Nevertheless it does not correctly account for observations of starburst- 
ing mergers: real mergers often have a nuclear starburst component, but also a 
relatively extended starburst activity, often taking place in massive star clusters. 
Recent simulations that resolve ISM turbulence and dense gas clouds show that 
the triggering of star formation in mergers is partly from gas inflows, but also 
from increased ISM turbulence, fragmentation into bigger/denser clouds, resulting 
in an excess of dense star forming gas compared to atomic gas or moderately-dense 
molecular gas. This can explain a "starburst regime" of star formation that does 
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Fig. 10. Kennicutt diagram comparing the surface density of gas and star formation in 
isolated disk galaxies (LTGs) and mergers near the peak of their starburst activity, as well 
as post-merger elliptical-like galaxies (ETGs). The "disks" and "starbursts" sequences 
are from Daddi et al. (2010). The low star formation efficiency in ellipticals is studied 
in l|Martig et all 120091 ) . 



not follow the global scaling relations observed for isolated disk galaxies. 



5 Structure formation in mergers, Super Star Clusters and Globular 
Clusters 

Galaxy interactions and mergers lead to bursts of star formation, but also play a 
major role in forming new stellar structures: rings, tails, shells, and more impor- 
tantly massive star clusters that could be the progenitors of globular clusters and, 
maybe in some cases, of tidal dwarf galaxies and/or ultracompact dwarf galaxies. 
Examples of formation of massive/compact star clusters in the tidal tails of major 
mergers are shown on Fig. [IT] A detailed review of structure formation and star 
cluster formation /evolution in galaxy interactions and mergers can be found in 
Bournaud! (|2010ali . 



6 Conclusion and implications for galaxy formation 

Galaxy interactions and mergers are a major driver of cosmic star formation, even 
if their contribution to the SF budget in the Universe compared to quiescent star 
formation in isolated disk galaxies is still unclear and debated. 

Gas compression in galaxy interactions and mergers results from two mecha- 
nisms: a global one (radial gas inflows towards the nucleus) and a local one (com- 
pression and fragmentation of the ISM in turbulent flows with a Mach number 
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well above one). The latter mode was really probed only recently, although obser- 
vations of triggered ISM turbulence in interactions have existed for more than a 
decade. Numerical simulations are hence barely starting to model merger-induced 
starbursts in a realistic way, and cosmological models probably do not account for 
these events in a realistic way, which will be a crucial aspect to understand the 
global star formation history of galaxies. 

Nevertheless, recent accurate observations of ISM properties and star formation 
in mergers can be accounted by theoretical models. These include a "double law" 
of star formation for mergers versus isolated disks, unveiling a "starburst mode" 
and a "quiescent mode", as well as an excess of dense molecular gas in the ISM 
of starbursting mergers. This probably shows that our theoretical understanding 
of the driving processes has improved and that modern simulations can model 
these in a realistic way. Major issues nevertheless remain, such as the role of 
stellar feedback processes, the triggering of gaseous outflows, the influence of gas 
accretion by central supermassive black holes in mergers, the mergers of such black 
holes into bigger and bigger black holes but with gravitational recoil effects, and 
the energetic feedback from these black holes. 
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